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ABSTRACT 

We present the discovery of two z > 6 quasars, selected as i band dropouts in the VST 
ATLAS survey. Our first quasar has redshift, z = 6.31 ± 0.03, z band magnitude, 
zab = 19.63 ± 0.08 and rest frame 1450A absolute magnitude, M 1450 = -27.8 ± 0.2, 
making it the joint second most luminous quasar known at z > 6. The second quasar 
has z = 6.02 ± 0.03, zab = 19.54 ± 0.08 and Mi 4 50 = -27.0 ± 0.1. We also recover a 
2 = 5.86 quasar discovered by Venemans et al. (2015, in prep.). To select our quasars 
we use a new 3D colour space, combining the ATLAS optical colours with mid-infra¬ 
red data from the Wide-field Infrared Survey Explorer (WISE). We use iab — zab 
colour to exclude main sequence stars, galaxies and lower redshift quasars, W1 - W2 
to exclude L dwarfs and zab — W2 to exclude T dwarfs. A restrictive set of colour 
cuts returns only our three high redshift quasars and no contaminants, albeit with a 
sample completeness of ~50%. We discuss how our 3D colour space can be used to 
reject the majority of contaminants from samples of bright 5.7 < z < 6.3 quasars, 
replacing follow-up near-infra-red photometry, whilst retaining high completeness. 

Key words: quasars: general - quasars: individual: ATLAS J025.6821-33.4627 - 
quasars: individual: ATLAS J029.9915-36.5658 


1 INTRODUCTION 

With the exception of transient events, bright quasars are 
the only objects luminous enough for high signal-to-noise 
ratio (SNR) spectra to be obtained out to very high red- 
shifts. Such observations are key to our understanding of 
both the evolution of the intergal actic medium (IG M) dur¬ 
ing the epoch of reionization (e.g. iFan et all 120061 ) and the 
growth of the first supermassive black holes. 

The Lyman alpha (Lya) forest observed in the spectra 
of quasars at moderate redshift (z « 3) arises from filaments 
of neutral hydrogen in the IGM. Models for reionization sug¬ 
gest that by z ~ 6, these filaments have started joining to- 
gether to form a neutral IGM w i th high filling factor (e.g. 
iKuhlen fc Faucher-Giguerell201 j . ISchroeder et alj|2013l ) re¬ 
sulting in the absorption of almost all flux at wavelengths 
blueward of Lya in the quasar rest frame and causing the 


* E-mail: a.c.carnall@durham.ac.uk 
f E-mail: tom.shanks@durham.ac.uk 


Gunn-Peterson troughs dGunn fc Petersonl 119651 ) observed 
in the spectra of quasars at these redshifts. 

Ionised regions at z > 6 are believed to be the direct 
result of photoionization by nearby AGN and star forming 
galaxies. Observations of these “bubbles” of ionised gas and 
the intervening neutral IGM, when illuminated by back¬ 
ground quasars, can be used to probe the metallicity and 
ionised fraction of gas in the IGM, and in the vicinity of pro- 
togalaxies during the epoch of reionization. In the case of the 
current high est redshift known qu asar, ULAS J1120+0641 
at z = 7.1 (IMortlock et alil201lh . the gas in the vicinity 
of the object was sho wn to be significant l y neutral and 
extre mely metal poor dSimcoe et al.l 120121 . iFinlator et al.l 
l2013h . providing a first glimpse of the formation processes 
influencing protogalaxies at z ~ 7, and of the evolution and 
enrichment of the IGM during reionization. 

Furthermore, the apparent existence of ~ 10 9 Mq black 
holes at this early point in cosmic history challenges our 
canonical pictures for the formation of black holes from ac¬ 
creting seeds, and the collapse of the dark matter haloes 
in which they reside. Two competing scenarios for mas- 
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sive bla ck hole format ion over such a short timescale ex¬ 
ist (e.g. H aima Firstly, the direct collapse of vast 

quantities of warm gas may take place, a scenario which 
would require very specific conditions for collapse in the 
early Universe. Secondly, the merger of many smaller black 
holes, formed from the first generation of population III 
stars may take place within deep potential wells due to 
massive dark matter haloes (> 10 12 Mq), potentially lead- 
ing to an anti-hierarchic al distribution of black hole masses 
dVolonteri fc Reesll2006l '). 

Future progress in these fields will require high res¬ 
olution spectroscopy of larger samples of z > 6 quasars. 
However, the spatial density of bright high re dshift quasars 
is ver y low, with the luminosity function of IWillott et all 
d20ldf predicting only ~110s>5.7 quasars with zab <20.0 
across the whole sky. 

Conveniently, the sharp drop in flux blueward of Lya 
due to IGM absorption allows high redshift quasar candi¬ 
dates to be selected from wide field optical survey data as 
dropouts at all wavelengths shorter than the z band for 
redshifts 5.7 < 2 < 6.5 and the Y band for objects with 
6.5 < 2 < 7.2. Howev er candidate lists obtained using this 
method, pioneered by |hjm_et_alJ d.2001 ) , are polluted by L 
and T dwarf stars dKirkpatrick et alT 1999l l which are far 
more numerous than high redshift quasars. 

This contamination means one or more further dimen¬ 
sions must be added to the colour space in order to dif¬ 
ferentiate between high redshift quasars and cool dwarfs. 
Initially follow-up J band photometry on a large num- 
ber of i band d ropout candidates was performed, e.g. by 
iFan et a hl d200l| and late r papers'). I Jiang et all (l2008l . [20081 1 
and Willott et aid d2007l . 1 20091 . l201ol l. These authors select 
objects with blue colours in zab — J to exclude L and T 
dwarfs. 

Where near-infra-red photometry is available from sur¬ 
vey data, a multidimensional colour spac e can be con¬ 
struct ed without follow-up observations, e.g. iBanados et aid 
d20l4) who use an i ab —zab vs zab —Y colou r space to select 
a popu l ation of 5.7 < 2 < 6.2 quasars, and lMortlock et aid 
(2009b|, 201l|) w ho us e a statis tical approach, d e scribe d in 
iMortlock et aid d2012li . Finally, IVenemans et al.1 (I2013T ) use 
a zab — Y vs Y — J colour space to select 2 > 6.5, z band 
dropout quasars. All of these methods however still require 
deep follow-up photometric observations in the i band to 
confidently exclude cool dwarf stars before proceeding to 
spectroscopic confirmation of their quasar candidates. 

In this work, we propose a novel method of three dimen¬ 
sional optical plus mid-infra-red colour selection to define 
an extremely clean sample of 5.7 < 2 < 6.3 quasars, using 
data from the Wid e-held Infrared Survey Explorer (WISE, 
iWriaht et aid l20ld i. detailed in Section El We therefore 
show that follow-up near-infra-red observations can be re¬ 
placed with publicly available mid-infra-red data to select 
the brightest i band dropout high redshift quasars in optical 
surveys. We then apply our method to the VST ATLAS sur¬ 
vey, within which we select two new 2 > 6 quasars directly 
from the survey data without the need for imaging follow-up 
observations. 

The structure of this paper is as follows. In Section[2] w e 
introduce the VST ATLAS survey, then set out the initial se¬ 
lection criteria we implemented to identify i band dropouts 
in our survey data. We then detail the optical plus mid- 
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Figure 1 . Colour diagram for separating dwarf stars from high 
redshift quasars. The shaded region denotes our selection criteria. 
L and T dwarfs are plotted in green and red respectively. The 
quasar trace based solely upon the SDSS composite is plotted in 
brown and the truncated, power law extended trace in black (see 
Section f2.3ll . Both have points at A 2 = 0.1 intervals, starting from 
z = 5.7. The 13 published quasars with z > 5.7, zab < 20.0 and 
Wl, W 2 SNR > 3.0 are plotted in blue. The three high redshift 
quasars detected in VST ATLAS (Tabic [T]) are plotted in pink. 


infra-red selection methods used to obtain a clean sample 
of zab < 20.0 high redshift quasars. In Section [3] we report 
the discovery of two bright quasars selected from VST AT¬ 
LAS survey data. In Section [5] we discuss the utility of our 
selection method and how it can be more widely applied to 
select i band dropout quasars from optical survey data. We 
conclude and summarise our results in Section 0 

Throughout this paper, the following cosmological pa¬ 
rameters are ass umed: Hp = 70 km s - 1 Mpc -1 , Qm = 0.28 
and 12 a = 0.72 dKomatsu et al.ll2009l h All magnitudes are 
quoted on the AB system except for WISE magnitudes, W 1 
(3.4 fim) and W 2 (4.6 fim) which are on the Vega system. 


2 CANDIDATE SELECTION 
2.1 The VST ATLAS Survey 

The Very Large Telescope Survey Telescope (VST) AT¬ 
LAS survey is an optical ugriz survey on the 2.6m VST 
at Paranal, which aims to image ~5000 square degrees of 
the Southern sky at similar depths to the Sloan Digital Sky 
Survey (SDSS) in the North. However VST ATLAS is deeper 
in the z band, with a mean 5er limiting magnitude of 20.89 
(as opposed to ~20.5 for SDSS), and also has considerably 
better seeing, wit h median seei ng in the range ~ 0.”8 — l."0 
for all five bands (IShanks et al.ll2015l h 

All VST ATLAS data is processed by the Cambridge 
Astronomical Survey Unit (CASU) and the calibrated im¬ 
ages and single band catalogues made available to the con¬ 
sortium and the Wide Field Astronomy Unit (WFAU) in 
Edinburgh. WFAU then produce band merged catalogues, 
available through the Omegacam Science Archive (OSA). 
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Quasar 

Redshift 

i (AB mag) 

z (AB mag) 

W 1 (Vega mag) 

W 2 (Vega mag) 

ATLAS J025.6821-33.4627 

6.31 ± 0.03 

>22.2 

19.63 ± 0.06 

16.12 ± 0.05 

15.48 ± 0.09 

ATLAS J029.9915-36.5658 

6.02 ± 0.03 

>22.7 

19.54 ± 0.08 

16.70 ± 0.08 

15.87 ± 0.13 

VIKINGKiDS J0328-3253 

5.86 ± 0.03 

>22.5 

19.75 ± 0.12 

16.90 ± 0.07 

16.19 ± 0.13 


Table 1. VST ATLAS and ALLWISE magnitudes for the three quasars selected from the VST ATLAS survey using the method described 
in Section [3] None of the objects are detected in the i band and 3a limiting magnitudes are provided. 


ATLAS DR1, a public release containing data up to 30th therefore decided to apply our selection to remove dwarf 

September 2012 was delivered to ESO in November 2013. stars before inspecting cutout images for the full sample. 


2.2 Initial Criteria 

For this paper we use the ATLASv20131029 proprietary re¬ 
lease, including data up to 31st March 2013 and containing 
2060 square degrees imaged in the i and z bands. To pro¬ 
duce an initial list of i band dropout objects within this area 
from the WFAU band merged catalogue, we specified that 
objects must have: 

(i) 18.0 < z AB < 20.0. 

(ii) Point source classification in the CASU catalogue 
(see Chehade et al. 2015, in prep, for a discussion of the 
robustness of this classification), and no WFAU error flags. 

(iii) No catalogued detection in the i band (the catalogue 
detection limit is 5<r) or iab — zab > 2.2. 

(iv) No catalogued detection in the u, g or r bands. 

The upper limit imposed by (i) is necessary in order to con¬ 
firm dropout status for objects undetected in the i band, the 
5cr i band detection limit for VST ATLAS being iab = 22.0. 
The lower limit on magnitude imposed by (i) w as based upon 
the pr ediction by the luminosity function of IWillott et al.l 
(120101) that the spatial density of z > 5.7, zab < 18.0 

quasars is less than one per 10,000 square degrees, therefore 
none were expected in our search area. Criterion (iv) ex¬ 
cludes a very small number of sources which are assumed to 
be either due to imperfections in the catalogue band merging 
process or poor i band data in some fields. 

After imposing these criteria on the data set, 9,607 ob¬ 
jects remained in our sample. In order to exclude any spuri¬ 
ous z band detections polluting the sample, and to facilitate 
our subsequent WISE colour selection process the following 
steps were then implemented. Objects must have: 

(v) A detection in the ALLWISE catalogue within 3.”0. 

(vi) SNR > 3.0 in W1 and W 2, corresponding to mean 
mag limits at -35° declination of 17.9 in W1 and 17.0 in W 2. 

After these criteria had been implemented, 3,452 objects re¬ 
mained in our sample. A visual inspection of a random sam¬ 
ple of these candidates revealed significant contamination by 
objects which were either misclassified as i band dropouts 
due to imperfections in the band merging process or which 
could not be confirmed as i band dropouts due to areas of 
poor i band data. 

To reject these objects, it is necessary to inspect cutout 
images of every candidate, however our sample is also still 
contaminated by L and T dwarfs. Our intention was to re¬ 
ject dwarfs using optical plus mid-infra-red colour selection 
and it was likely that a significant proportion of the other 
contaminants would also be rejected by this process. It was 


2.3 Optical Plus mid-infra-red Selection 

It has been well documented that quasars have redder colour 
in W1 - W 2 than main sequence stars fe.g. lStern et al.l|2012l . 
iBlain et al. IHoH). L dwarf stars also have bluer W1 - W2 
colours, which we can exploit to reject L dwarfs from i band 
dropout samples by a W1 - W2 colour cut. T dwarfs in¬ 
habit a temperature range for which the black body com¬ 
ponents of their spectra peak in the mid-infra-red, and also 
exhibit temperature dependent methane absorption in the 
W1 band, giving them a broad range of redder W1 - W2 
colours. However these objects are very dim by comparison 
in the optical, allowing us to split them off from high redshift 
quasars up to z ~ 6.3 by their very red zab — W2 colour. 

In order to define suitable colour cuts, we first collected 
catalogues of L and T dwarfs in order to assess the dis¬ 


stars from Burgasser et al. (2004). Scholz et al. (20091). 

Schmidt et al. (2010]). Kirk 

patrick et al. (2010). Revie et al. 

hold). Albert et al.l (2011 

) and Kirkpatrick et al.l 201 ll) 


and matched them with the ALLWISE catalogue. 

To understand the distribution of high redshift quasars 
in colour space we obtained magnitudes for all 17 pub¬ 
lished quasars with z > 5.7, zab < 2 0. 0 jFan et al.ll2006l. 
Mortlock et alJ l2009bl . iBanados et al.l 12014 and Wu et all 


20151) . Of these. 13 have SNR greater than 3.0 in the W1 


and W2 bands. We also produced a trace of the evolu¬ 
tion of quasar colour with redshft by applying synthetic 
photometry to the media n composite quasar spectrum of 
IVanden Berk et al.l d200lTl (hereafter the SDSS composite). 

IGM Lya absorpt ion was modelle d by applying the cor¬ 
rections proposed by ISongalid (120041 ) between rest frame 
1026A and 1216A. Flux blueward of Ly/3 was assumed to be 
zero. The SDSS composite is known to be subject to greater 
host galaxy contamination at rest frame wavelengths above 
~5000A than is expected for highly luminous high redshift 
quasars. We attempted to mitigate this effect by truncating 
the SDSS composite at rest frame 4200A and extending to 
longer wavelengths using the bluer a v = -0.44 power law 
used to fit the SDSS composite bluewards of 5000A. Figure 
|T] shows the W1 — W2 vs zab — W2 plane with the resulting 
trace plotted in black, along with the other data samples. A 
trace produced without truncation of the SDSS composite 
is shown in brown for comparison. 

For this work, we defined a restrictive set of colour cri¬ 
teria based upon this information to select the best quasar 
candidates in our field. Our final criteria are as follows. Ob¬ 
jects must have: 


(vii) 0.55 < W1-W2 < 1.2 
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Figure 2. LDSS-3 spectra for the two newly di scovered quasars. The spectra have been binned to a resolution of 2 A. Data reduction was 
carried out using the Astropy Python package <|Astropy Collaboration! 12013l ). Instrumental response was corrected using an observation 
of the standard star LTT 3864. Absolute flux calibration was carried out using synthetic photometry to scale to the z band magnitudes 
obtained from VST ATLAS. Positions consistent with our calculated redshifts are shown for the Ly-a, Ly-/3 and NV emission features. 


(viii) 3.2 < zab - W 2 < 4.5 

These colour criteria are also shown in Figurc[l] It was found 
that W1 - W2 colour greater than 0.55 is a good criterion for 
excluding L dwarf stars, and zab - W2 < 4.5 discriminates 
well between T dwarfs and high redshift quasars up to 2 ~ 
6.3. The upper limit on Wl - W2 and the lower limit on zab 
- W2 were chosen so as to reject objects misclassified as i 
band dropouts and thus reduce the number of candidates to 
be visualy checked, without rejecting objects with colours 
similar to high redshift quasars. It must be stressed that 
the colour selection shown in Figure |T] is only effective at 
selecting high redshift quasars if the sample being considered 
are all i band dropouts. The selection of Figure [I] alone is 
not enough to exclude extragalactic contaminants such as 
lower redshift quasars. 

After applying (vii) and (viii), 200 objects remained 
within our selection criteria. Upon visual inspection using 
the CASU postage stamp browser, 197 of these were found to 
be either obvious misclassifications for which the i band de¬ 
tection had not been successfully associated with the z band 
detection, or areas of poor i band data where dropout status 
could not be confirmed. Magnitudes for the three remain¬ 
ing sources are given in Table [I] Of these sources, one was 
recognised as the only known z > 5.7, zab < 20.0 quasar in 
our search area, VIKINGKiDS J0328-3253 (Venemans et al. 
2015, in prep.). As the luminosity function of lWillott et al.l 
d 2010 T) predicts six z > 5.7, zab < 20.0 quasars within our 
search area, we decided to proceed to spectroscopic confir¬ 
mation of the two remaining candidates. 

3 TWO NEW Z > 6 QUASARS 

3.1 Spectral Observations 

We obtained long-slit spectra for our two candidates on the 
night of 16th January 2015 using the LDSS-3 instrument 


on the Magellan-II telescope at Las Campanas Observatory 
(observer M. Rauch). Exposures of 600s for each object were 
taken using the VPH RED grism and a l.”0 slit. The data 
obtained are shown in Figure [5] Both objects show a strong 
spectral peak in the near-infra-red which we identify as Lya 
emission, and a break in continuum flux blueward of this 
point. We now discuss the features of each spectrum in turn. 


3.2 ATLAS J025.6821-33.4627 

This quasar exhibits a broad Lya emission line from which 
we have estimated a redshift, 2 = 6.31 ± 0.03. By scaling the 
template spectrum employed in Section \2. 31 to the observed 
redshift and z band magnitude of the object we estimate 
it has M 1450 = -27.8 ± 0.2. This is the s ame as estimated 
for SDSS J1148+5251 bv lFan et al] J2006lb making this the 
joint second most luminous k nown z > 6 quasar, after SDSS 
J0100+2802 rtWu et al.ll2015h . 

The spectrum exhibits a faint signal consistent with 
Ly /3, and an excess in transmission between 8700A and 
8800A. This could indicate an ionised region close to the 
quasar, however near-infra-red spectroscopy is needed to 
better constrain the exact quasar redshift and hence the 
nature of any ionised regions. 


3.3 ATLAS J029.9915-36.5658 

For this quasar, we estimate a redshift of 2 = 6.02 ± 0.03 
based upon the position of the broad Lya line. Following the 
same procedure as above, we estimate an absolute magni¬ 
tude of M 1450 = -27.0 ± 0.1. There is a probable absorption 
feature superimposed on the Lya emission line, as can also 
be observe d in the spectra of se veral of the quasars first 
reported bv lBanados et aD d2014l 'l. 
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4 DISCUSSION 

In this work we have employed a restrictive set of colour 
criteria. This approach was chosen in order to prove the 
concept of our optical plus mid-infra-red selection method 
by reducing contamination to a very low level. Naturally 
this takes a toll on the completeness of the sample of 
high redshif t quasar s obtained. The luminosity function of 
IWillott et al.l (l2010l l predicts six quasars with redshift 0 

> 5.7 and z band magnitude zab < 20.0 over the 2060 
square degrees of our search area. The confirmation of all 
three of our candidates as high redshift quasars confirms the 
extremely clean nature of our colour selection, and places 
the completeness of our sample at ~50%. This is supported 
by the fact that 9 of the 17 published 2 > 5.7, zab < 20.0 
quasars fall into our selection criteria, suggesting a complete¬ 
ness of ~53%. 

When considering how to extend our colour selection to 
obtain a more complete sample of quasars, it is instructive 
to analyse the location in colour space of the eight quasars 
which do not meet our selection criteria as plotted in Figure 
[T] Four quasars are not detected in the ALLWISE source 
catalogue within 3.”0 with SNR above 3.0 in the W1 and 
W 2 bands. A further two fall to the left of our selection 
box, with W1 - W 2 colours of 0.38 and 0.44, closer to the 
distribution of L dwarf stars. Two other objects fall just 
above our selection box with zab — W2 = 4.56 and 4.69. 

Our selection criteria could be extended and coupled 
with follow-up observations in the optical to select a more 
complete population of bright 5.7 < 2 < 6.3 quasars. A nat¬ 
ural method for this would involve relaxing the lower bound 
on 171 - W2 and upper bound on zab — W2, then con¬ 
ducting follow-up i band observations to exclude L and T 
dwarfs within the extended selection region. This selection 
still results in the removal of the majority of L and T dwarfs, 
meaning a significant reduction in the number of candidate 
objects for which follow-up observations are required. 

Our method has two limitations, the first of which is the 
requirement that candidates be bright enough in the mid- 
infra-red to appear in the W 1 and W2 bands with a SNR 
greater than 3.0. Currently, this limits the usefulness of our 
method in selecting lower luminosity quasars, however on¬ 
going and future deeper surveys in the mid-infra-red such 
as the Spitzer-IRAC Equatorial Survey (SpIES) and NEO¬ 
WISE will soon allow this method to be applied to objects 
with lower intrinsic luminosities. Secondly the steep increase 
of the zab — IE2 colour of quasars with redshift above 2 ~ 
6.2 means we cannot select clean samples of quasars with 2 

> 6.3 using our method. An analogous colour selection for 
higher redshift objects using Y—W2 colour fails to efficiently 
separate T dwarfs from high redshift quasars. 


5 CONCLUSION 

In summary we have proven the concept of a new method of 
optical plus mid-infra-red colour selection for 5.7 < 2 < 6.3 
quasars. Using a restrictive set of colour criteria we identified 
two candidates which were confirmed to be bright quasars 
at redshifts of 6.31 ± 0.03 and 6.02 ± 0.03, with M 1450 = 
-27.8 ± 0.2 and -27.0 ±0.1 respectively. The former is the 
joint second most luminous known 2 > 6 quasar and both 


are well positioned for observation by powerful current and 
future Southern observatories. Our method allows selection 
of i band dropout quasars to proceed from optical catalogues 
with a significant reduction to the volume of follow-up ob¬ 
servations by replacing Y or J band follow-up photometry 
with publicly available mid-infra-red WISE data. 
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